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As part of an international Rn mea- 
surement intercomparison conducted at 
Bermuda in October 1991, NIST provided 
standardized sample additions of known, 
but undisclosed ("blind") ^^^Rn concen- 
trations that could be related to U.S. 
national standards. The standardized sample 
additions were obtained with a calibrated 
^^^Ra source and a specially-designed 
manifold used to obtain well-known dilu- 
tion factors from simultaneous flow-rate 
measurements. The additions were intro- 
duced over sampling periods of several 
hours (typically 4 h) into a common 
streamline on a sampling tower used by 
the participating laboratories for their 
measurements. The standardized ^^^Rn 



activity concentrations for the intercom- 
parison ranged from approximately 
2.5 Bq • m"'' to 35 Bq • m"' (of which 
the lower end of this range approached con- 
centration levels for ambient Bermudian 
air) and had overall uncertainties, approxi- 
mating a 3 standard deviation uncertainty 
interval, of about 6 % to 13 %. This paper 
describes the calibration and methodol- 
ogy for the standardized sample additions. 
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1. Introduction 

In October 1991, an international measurement inter- 
comparison of instruments used to measure trace atmo- 
spheric concentrations of ^^^Rn (radon) was organized 
by Drexel University and conducted at the Atmosphere/ 
Ocean Chemistry Experiment (AEROCE) test site 



located at Tudor Hill, Bermuda which is operated in 
conjunction with the Bermuda Biological Station for 
Research, Inc. This intercomparison comprised two 
parts: Viz. (1) measurement comparisons among four 
laboratories of commonly sampled ambient air over 
approximately a 2 week test period; and, (2) measure- 
ment comparisons among three of these four laborato- 
ries of a select number of introduced samples with 
known radon activity concentration. 

Measurements of atmospheric radon in remote 
marine environments are used to obtain information on 
the temporal and spatial distributions, which are in turn 
used to test and validate global models that simulate the 
transport and removal of trace atmospheric species from 
continental air masses. Unlike other chemical species, 
radon is an excellent tracer for such studies because it 
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has a well-characterized source (large land masses from 
which radon diffuses ubiquitously) and only one princi- 
pal "sink" (radioactive decay). For such measurement 
data to be useful, however, particularly when collected 
by different laboratories at different global sites using 
diverse instruments that are based on widely different 
collection and measurement principles, it is necessary 
that the collected data have at least a common relative, 
if not an absolute, reference basis. The intent of this 
intercomparison was to provide that basis for the partic- 
ipating laboratories. 

The participating laboratories in the intercomparison 
were: Environmental Measurements Laboratory, U.S. 
Department of Energy (New York, USA); Drexel Uni- 
versity (Philadelphia, USA); Australian Nuclear Sci- 
ence and Technology Organization (Menai NSW, Aus- 
tralia); and the Centre des Faibles Radioactivitis, 
Laboratoire Mixte C.N.R.S. -C.E.A. (Gif-sur-Yvette, 
France). Henceforth referred to as Lab E, Lab D, Lab A, 
and Lab F, respectively. The first three laboratories per- 
formed simultaneous measurements from a common 
streamline on an ambient air sampling tower. Lab F, in 
contradistinction, sampled ambient air nearly adjacent to 
the inlet of the sampling tower and therefore only partic- 
ipated in the ambient air intercomparisons. 

The National Institute of Standards and Technology 
(NIST) was invited to participate in this intercompari- 
son to serve as the ' 'referee' ' and, if possible, to provide 
an "absolute" (sic) reference basis for the participants' 
measurements. The planned scheme for the latter role 
consisted of NIST being able to provide standardized 
sample additions of known, but undisclosed ("blind") 
radon concentrations that could be related to U.S. na- 
tional standards. The standardized sample additions 
were to be introduced over sampling periods of several 
hours into a common streamline on the sampling tower 
used by the participants for their measurements. Beside 
being a "blind" intercomparison in terms of the radon 
activity concentrations, the timing and duration of the 
standard sample additions were also largely unknown to 
the participating laboratories. Radon concentrations for 
the sample additions were planned to vary from approx- 
imately a few times typical ambient concentration levels 
up to about 1000 times ambient. This necessitated being 
able to provide well-known sample dilutions over a wide 
dynamic range. 

The sample additions were made with a commer- 
cially-available, solid ^^''Ra source, calibrated by NIST 
in terms of the available ^^^Rn concentration as a func- 
tion of flow rate through the source. The source was 
employed in conjunction with a sampling and dilution 
gas-handling manifold also provided by NIST. "Grab 



samples' ' from the sampling streamline were also taken 
with sample bulbs that were returned to NIST for assay- 
ing the ^^^Rn concentration. These "grab samples" 
were taken as confirmatory measurements to insure that 
the source and manifold were not performing differently 
at the test site as at NIST where the source calibrations 
were made. 

This paper summarizes and describes in detail the 
experimental arrangements, the NIST calibrations, 
methodology for the standardized additions, and the pro- 
tocol for the experimental aspects of the intercompari- 
son. The results of the intercomparison measurements 
for the participating laboratories are presented in a sep- 
arate companion paper [1]. 

2. Instrumentation and Methodology 

The experimental arrangement used for the intercom- 
parison is illustrated in the schematic diagram. Fig. 1. 

Ambient air was sampled continuously from the top 
of a sampling tower located at the Tudor Hill, Bermuda 
test site. The sampling stack consisted of a nominal 10.2 
cm diameter aluminum tube whose inlet was approxi- 
mately 20 m above ground level. The inlet at the top of 
the tube was a goose-neck, which was inverted to min- 
imize the intake of rain, and covered with a plastic mesh 
to prohibit the intake of birds, insects, etc. A similar 
diameter polyvinyl chloride (PVC) tube was coupled to 
the aluminum tube about 2.5 m above a horizontal sec- 
tion of the stack that was used by the participants. Air 
flow within the tube was maintained by sampling pumps 
that were part of the measurement instrumentation for 
three of the participating laboratories (Labs A, D, and 
E). These instruments were located at the terminus of 
the horizontal section. One laboratory (Lab F) per- 
formed measurements at the top of the sampling tower 
nearly adjacent to the inlet and, as mentioned previously, 
did not participate in the intercomparison of standard- 
ized additions. The standardized additions of known 
^^^Rn activity concentration, and with known dilution 
factors, were made into the sampling stack at ground 
level at the foot of the sampling tower as shown in 
Fig. 1. 

The NIST manifold was designed to obtain known 
radon concentrations from a calibrated ^^''Ra source, to 
provide standardized dilutions of this concentration, and 
to introduce them into the sampling streamline. The 
manifold also allowed the filling of sample bulbs that 
were used for confirmatory measurements of the pro- 
vided radon concentration. The manifold was con- 
structed, for the most part, of 1/4 in and 3/8 in 
(0.635 cm and 0.953 cm) stainless steel (SS) tubing with 
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Fig. 1. Schematic diagram of the experimental configuration used for the measurement intercomparison 
showing the NIST manifold that was employed to provide standardized additions of known radon concen- 
tration and the relative sampling locations for the participating laboratories. 



compression Swaglok' and IPT threaded pipe fittings. 
All of the valves were of high-vacuum, stainless steel 
construction. The two pumps were simple, vibrating-di- 
aphragm, aquarium pumps which had maximum flow 



Certain commercial equipment, instruments or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the purpose. 



rates of approximately 1.5 L- min"^ and which utilized 
ambient air outside their housings. The two control 
valves (VC in Fig. 1) located after the source pump were 
used to vary and control the flow rate through the 
radium source. With the use of appropriate triad of 
valves, the flow rate could be adjusted with the stream- 
line bypassing the source. A similar triad of valves was 
used at the sampling port to either bypass the port or to 
direct the streamline of air through a sample bulb. 
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The ^^''Ra source was a commercially-available, flow- 
through source which was independently calibrated by 
NIST and will be subsequently described in Sec. 3 of 
this paper. The flow rates through the radium source (f4) 
was measured with flow meter F4. The streamline after 
meter F4 could be either bypassed to a waste stream or 
directed to the sampling stack for a standard addition. 
The flow rate (f^) from the dilution pump was adjusted 
with a single, control valve teed off its streamline. This 
rate was measured with meter F3. Flows /t and/3 were 
added and mixed through a tortuous mixing path that 
consisted of a crossed path of flexed SS bellows tubing 
which served as a flow turbulator. Its efficacy and the 
adequacy of mixing were demonstrated by confirma- 
tory measurements that will be reported in Sec. 5 of this 
paper. A portion of the combined (/I, + fj) flow rate 
could be diverted through a control valve to a second 
waste line. The net flow rate (f^) after this diversion was 
measured with meter F5. The three meters (F4, F3, and 
F5) where matched and intercalibrated mass flow 
meters (Matheson Gas Products, East Rutherford, NJ) 
having a range of L • min"' to 3 L • min"' with 
standard volume conditions of 21 °C and 760 mm 
Hg (1.013 X 10' Pa), and having a reported calibration 
accuracy of 1.0 %. Confirmatory measurements of their 
in situ intercalibration are also reported in Sec. 5 of this 
paper. The two waste streamlines, which at times con- 
tained appreciable radon concentrations, were directed 
by flexible, plastic hose to distances considerably away 
from the sampling stack. For the first few standard addi- 
tions, the waste hose outlet was located 30 m to 40 m 
southeast of the tower base, and was subsequently re- 
located about 30 m down the hill from the tower (which 
would place the outlet 40 m to 50 m below and west of 
the intake at the top of the tower). None of the diversions 
were considered to have appreciably altered any ambient 
radon concentration levels. 

The outlet of the NIST manifold with flow rate/5 was 
directed into the main sampling stack where the flow of 
the combined streamlines was measured with meter F 1 
and mixed with a commercial, high-efficiency flow 
tubulator prior to sampling by the instrumentation of 
Labs A, D, and E. As specified by the manufacturer, 
meter Fi was a ft^ • min"' to 50ft^-min"' 
(0 L • min"' to 1416 L • min"' mass flow meter (TSI, 
Inc., St. Paul, MN) with standard volume conditions of 
70 °F (20.1 °C) and 14.7 lb • in"' (1.013 X 10' Pa) and 
reported "tolerance" of ±2 % of reading plus 0.2 % of 
full scale. 

Throughout the intercomparison of standardized 
additions, all four flow meters were continuously 
monitored. Their V to 5 V analog outputs were 



sampled at 10 s intervals and converted to digital signals 
with a 4096 channel analog-to-digital converter (ADC) 
and chronologically recorded on a dedicated personal 
computer (PC). This every 10 second file of all of the 
flow rates provided an exhaustively complete record of 
all flow-rate changes, and thereby the necessary dilution 
factors, throughout the course of all the standardized 
additions. The chronological record was maintained 
with the internal computer clock. This clock was curso- 
rily intercompared with the timing devices of the 
measurement instruments of the participating labora- 
tories. All clocks were in agreement within a minute or 
two, and any of the differences in clock times could be 
demonstrated to have negligible effect on the intercom- 
parison results. 

A simplified schema of the flow network is shown in 
Fig. 2. The radon concentration in the source and 
dilution lines are 



and 



C4 — Cs + C . 



C3 — Ca : 



(1) 



(2) 



respectively, where Ca is the '^'Rn activity concentra- 
tion in ambient air and Cs is the '^'Rn activity concentra- 
tion from the radium source when maintained at flow 
rate/4 (see Sec. 3). With C4 at flow rate/4 and C3 at flow 
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Fig. 2. Simplified flow network for the experimental configuration 
shown in Fig. 1. 
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rate f^, the concentration at their confluence is just the 
weighted sum: 



C 



J 4 C 4 + /3 C 3 
(A+/3) 



(3) 



/4 



/4+/3 



Cs + Ca . 



(3) 



With a partial diversion of the combined flow rate 
(f4 + /s) to /x and fs, the concentration exiting from the 
NIST manifold streamline is diluted to 



C, = 






(1 



c + 



/.-/5 



Ca 



(4) 



/4+/3 



Cs+Ca 



The protocol for initiating the standardized additions 
consisted of: (1) adjusting flow-rate /t while bypassing 
the radium source; (2) opening the flow /4 through the 
source but diverting it out the waste line to allow time to 
reach a steady-state concentration Cs from the source 
and to flush the source line; (3) adjusting flow rates /s 
and/5 with the appropriate control valves; and finally, 
(4) after all the flow rates are stabilized, and at the 
commencement of the standard addition, opening the 
valve after meter F4 and closing the waste bypass valve 
to direct Cs into the manifold stream. The standardized 
additions were typically of 4 h duration (with one excep- 
tion) to allow sufficient overlap with the measurement 
intervals for each of the participating laboratories. Each 
addition was initiated at the start of one of the hourly 
measurement cycles of Lab E. The rapid and abrupt 
change in flow-rate /i at those hourly intervals (a flow- 
rate change of approximately 400 L ■ min"' over approx- 
imately 40 s as the Lab E pump went off and on again) 
provided a convenient internal timing mechanism. The 
standard addition was terminated, again at the end of 
one of the measurement cycles of Lab E, by : ( 1 ) rapidly 
diverting Cs from the main streamline to the bypass 
waste line; and (2) readjusting the flow rates (and 
thereby the dilution factor) for the next standard 
addition. 

3. Calibration of the ^^"^Ra Source 



in the main sampling line. This is the radon concentra- 
tion sampled by the participating laboratories. For dis- 
cussion purposes, let us call the portion contributed by 
the standard sample addition Co, and call the combined 
flow-rate ratios the dilution factor D , so that: 



For a ^^^Ra source having activity Ar^, the instanta- 
neous rate of production of ^^^Rn activity, Ar„, is simply 
given by 



dA. 



dt 



■ = R = Ar, All 



(8) 



and 



C] — Co + Ca ; 
Co = DCs ; 



D 



h 



/4+/3 



(5) where Ar„ = 2.09 X 10"^ s"' is the decay constant for 
^^^Rn. This result may be derived from several straight- 

(6) forward approaches, including taking the time-differen- 
tial limit of the generalized Bateman equation steady- 
state solution for the growth of a radioactive daughter 
from decay of its parent. 

Now, for an idealized flow-through ^^^Ra source, 
under continuous constant-flow conditions with com- 

(7) pletely mixed, turbulent air, the expected ^^^Rn activity 
concentration Cr„ in the source's exiting streamline 
would be 



It should be noted that there are two simplified situa- 
tions: with/t = and/3 = 0,/4 =fi and the dilution factor 
is just D = (fs/fi) = (/4//1). Alternatively, with /, = 0, 
(/s +/4) =/5 and D = (/4//1). For confirmatory 
purposes, both cases were tested in the series of stan- 
dardized additions (see Sec. 5). 



_ 0_ 

Crp — JV „— Ar^Arh „ 



(9) 



where (p is the emanation fraction (i.e., the fraction of 
the total ^^^Rn generated in the source that is released 
into the air stream) and / is the flow rate through the 
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source. This approach of obtaining the radon activity 
concentration in the streamline by dividing the contin- 
uously available radon activity by a corresponding flow 
rate is, of course, simplistic. Firstly, it ignores any time 
dependence necessary to reach a steady-state concentra- 
tion Cr„ in the streamline. This, however, is not overly 
problematic since a steady-state concentration is 
reached within a few minutes of continuous flow in 
most situations involving a reasonable range of flow 
rates and streamline volumes. It merely requires that 
sufficient time elapse to remove by flow ventilation any 
previously accumulated radon when the source was 
closed or operating at a different flow rate. Secondly, the 
approach of Eq. (9) assumes that the rate of removal of 
radon from the source by flow ventilation (or dilution) 
is rapid compared to its rate of generation. The rate 
constant A, for the rate of removal by ventilation in 
completely mixed, turbulent air is the quotient of the 
flow rate /by the internal volume v of the source, i.e., 
Av =//v. Again, for any reasonable range of flow rates 
and internal source volumes A, » Ar„, and the condi- 
tional is easily satisfied. (It may be useful to note in 
passing, however, that for users of ^^^Th/^^'^Ra sources 
to produce short-lived ^^"Rn ("thoron") with decay 
constant A = 0.0125 s"', the rate constants A and Av are 
often of comparable magnitude. In this case, all of the 
rate-of -change Bateman balance equations used to 
derive the steady-state solutions need to be modified by 
substituting (A + A,) for A.) 

The simple model of Eq. (9) implies that the steady- 
state ^^^Rn concentration in the streamline, Cr„, is just 
inversely proportional to the flow rate / through the 
source. The model, of course, is valid only if the emana- 
tion fraction <p is independent of flow rate, as well as of 
any other external variables. In addition, one can intu- 
itively imagine that there must be only some reasonable 
operating range of flow rates for any realistic source. 
One could conceive that at exceedingly low flow rates, 
not all of the radon at the emanating surfaces of the 
source would be adequately entrained into the stream- 
line for transport out of the source at constant time rates. 
Thus, because of possible variations in <p with flow rate 
or other external variables, as well as possible losses due 
to non-uniform transport at low flow rates, a calibration 
curve consisting of Cr„ as a function of flow rate/might 
very well exhibit deviations away from a strict 1 // pro- 
portionality. 

The flow-through radium source used for the inter- 
comparison was calibrated in terms of the available 
^^^Rn concentration Cs maintained in a streamline at a 
constant flow rate /4. For these calibration measure- 
ments, the streamline was sampled by filling flow- 
through sample bulbs. The source was operated and the 
streamline sampled using the upper portion of the 



identical manifold described above in Sec. 2. This 
portion consisted of the source pump, appropriate 
control valves to adjust the flow rate through the source, 
a triad of valves to serve the inlet and outlet ports for the 
source and a bypass line, another triad of valves for the 
sampling ports, and flow meter F4. 

The sample bulbs which were filled at known flow 
rates /t were assayed for total ^^^Rn activity using the 
NIST primary ^^^Rn measurement system. This system 
is based on pulse-ionization-chamber (PIC) measure- 
ments calibrated against national ^^^Ra solution stan- 
dards maintained by NIST, and has been described in 
detail previously by Colle et al. [2] and Hutchinson et al. 
[3]. The estimated uncertainty in the measurement of 
the total activity in the sample bulbs (corresponding to 
an approximate relative standard deviation) is approxi- 
mately 0.7 % [3]. 

The volumes of the sample bulbs, used to obtain the 
activity concentration Cs from the total activity in the 
samples, were well known from replicate gravimetric 
determinations using high-purity water of known den- 
sity at given temperatures. The sample bulbs were of 
various sizes ranging in volume from about 330 mL to 
670 mL. The estimated overall relative standard uncer- 
tainty in the volumes is 0.2 % to 0.5 %, the magnitude 
of which was primarily controlled by the number of 
replicates. 

Samples were collected and subsequently assayed for 
^^^Rn to determine Cs over a range of flow rates from 
about 0.3 L • min"' to 1 .3 L • min"'. Continuous flow for 
a minimum of 10 min to 20 min (i.e., at least 10 com- 
plete air changes through the manifold) was maintained 
during sampling to flush the manifold and remove all 
previous air, and to insure adequate mixing. The preci- 
sion of the flow-rate measurements during sampling 
was typically 0.02 % (relative standard deviation of the 
mean). Flow meter F4 was calibrated by its manufac- 
turer (Matheson) and was stated to have a calibration 
accuracy of 1 .0 % for air over the entire L • min"' to 
3 L • min"' flow rate range. Inasmuch as this uncer- 
tainty is presumed to correspond to some type of 
"maximum" confidence interval, the approximate 
relative standard deviation was taken to be 0.5 %. 

The calibration results for 16 independent determina- 
tions of Cs as a function of mean/4 are shown in Fig. 3. 
The data was fitted to obtain the ;^f ^-minimized regres- 
sion curve 



Cs 



: a + 



/4 



(10) 



with 



-2.386 (1 ± 16.0 %) Bq • L" 



Volume 101, Number 1, January-February 1996 

Journal of Research of the National Institute of Standards and Technology 



50 
45 
40 
35 
30 
25 
20 
15 
10 
5 



; r-' 1 1 r-n 1 , ^ ^-. . 

' *♦ '. 

L ♦ J 

r ♦ 

" ♦■ . ": 

^ "■♦■♦■♦. ^ ^ 

r ♦■■■■«■ ♦ - 



Q.4 0.8 O.a 1.0 

Flow Rate (L/min) 



1.2 1.4 



Fig. 3. Calibration curve for the Ra flow-through source terms of the concentration of 
^^^Rn in a streamline, Cs, as a function of the flow rate in/4. 



and 



b= 14.825(1 ± 1.5%)Bq-L- 



when/4 is expressed in units of L • min"'. The fitted 
residuals ranged from 0. 1 % to 1 .4 % with a correlation 
coefficient of r = 0.9968. Regressions using other fit- 
ting functions, including second- and higher-order poly- 
nomials in I//4 as well as exponentials, were inferior to 
that of Eq. (10). The uncertainties in Cs (for the flow 
rate range /4 = 0.3 L • min"' to 1.3 L • min"') are esti- 
mated to be 1.2 % to 2.6 % (relative standard uncer- 
tainty). This was obtained by propagating in quadrature 
the combined standard uncertainty in the sample assays 
for ^^^Rn concentration (0.7 % to 0.9 %), the flow meter 
standard uncertainty (0.5 %), and the standard uncer- 
tainty in Cs from the regression (0.77 % to 2.4 %). 
Refer to Table 7 in Sec. 6 for further detail. 

The ^^"^Ra source was a "Type RN-1025 Flow- 
Through Radon Gas Source" manufactured by Pylon 
Electronic Development Company (Ottawa, Canada). It 
was certified by Pylon in 1981 to contain 115.9 kBq 
^^''Ra (decay corrected to October 1991) with a total 
calibration accuracy of ±4 % at the 99 % confidence 
level. The Pylon certificate also gave a continuously 
available ^^^Rn activity rate of 14.59 Bq • min"'; and 
indicated that in continuous flow the ^^^Rn activity 
concentration (in units of Bq • L"') from the source 
could be calculated using the formula 



C = 



14.59 



(11) 



over the flow rate (F) range of a few tenths of a 
L • min"' up to 10 L • min"'. Comparison of Eqs. (10) 



and (11) indicate that there can be substantial differ- 
ences between the calibration curve and that given 
by Pylon. Over even the small flow-rate range of 
0.3 L • min"' to 1.3 L • min"', the Pylon valve exceeds 
that from the calibration curve by 3.4 % at /4 = 0.3 
L • min"', by 1 1 % at/4 = 0.7 L • min"', by 17 % at/4 = 
1.0 L • min"', and by over 24 % at/4 = 1.3 L • min"'. 
These systematic differences have been consistently 
borne out by observations by NIST with this radium 
source over the past 8 years over even larger flow-rate 
intervals (up to 3.5 L • min"'). 

It must be emphasized, however, that the calibration 
function of Eq. (10) is applicable for only this particular 
source and not meant to be applicable to any others. The 
model for the calibration function is phenomenalisti- 
cally, not theoretically, based; it is meant to be practical 
for its intended use over the small flow-rate range. 
Nevertheless, the previous discussion suggested why 
there could be variations from the 1 //functionality. Even 
the functional form of Eq. (10) for this particular source 
is probably not adequate over a wider flow-rate range. 



4. Radon-222 Standardized Additions 

Fifteen sets of standardized additions using the 
methodology described in Sec. 2 were performed over a 
period of 5 d (October 9-13, 1991). The duration of each 
standard addition was 4 h, with the exception of # 1 3 
which was 3 h. The ^^^Rn activity concentrations, Co of 
Eq. (6), for the standard additions ranged nominally 
from 2.5 Bq • m"^ to 35 Bq • m"^ (refer to Table 1). One 
of them (#4) had to be discarded from the analysis 
because a bypass valve on the manifold was inadver- 
tently opened during the course of the addition which 
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diverted the source flow stream from the main sampling 
line. 

The standardized additions were analyzed in terms of 
calculating a mean Co individually for each partici- 
pating laboratory for each of their sampling/measure- 
ment intervals. For Lab E, the start and stop times of 
their 1 h sampling intervals were easily obtained from 
the abrupt changes in the flow rate/i data (see Sec. 2). 
For Lab A, whose instrument recorded continuously and 
gave averaged results that were "smoothed" by a time 
constant of approximately 90 min [1], the entire 4 h 
(or 3 h in the case of #13) standard addition time 
interval was used to derive the mean Co. For Lab D, the 
mean Co was calculated using the start and stop clock 
times for their 2 h sampling/measurement intervals 
which were submitted with their measurement results. 
As indicated earlier (Sec. 2), the flow-data clock time 
was cursorily intercompared and corrected to the 
standard (GMT) clock time used by Lab D. 

To obtain a mean Co for each standardized addition 
interval from some start time t^ to stop time f^ for each 
laboratories' individual measurement cycle, the flow 
data could be analyzed in at least three different ways but 
all yield essentially the same numerical results. 

The first approach consisted of calculating individual 
mean flow rates f 4, f^, /s, and /i, over the given 
fa to fz interval. The mean /4 was then used to obtain 
a mean Cs from Eq. (10); and a mean dilution factor 
D was obtained from the individual flow-rate means 
using Eq. (7). The mean Co was thereby approximated 
by the product of the derived Cs and D [Eq. (6)]. 



The second approach consisted of calculating the 
quantities Cs(ti), D(ti), and Co(f,) from the flow-rate 
measurement data/4(f, ),f?,{ti ),fs(ti ), and/i (f, ) at each 
ti measurement point (taken 10 s apart). The set of 
resulting Co(f, ) values was then arithmetically averaged 
over the entire t^ to t^ interval to obtain the mean Co. 

The third approach, considering the temporal varia- 
tion in Co, consisted of numerically integrating the 
individual Co(f,) values at each point f, normalized by 
the total flow (from flow rate /i ) in the main sampling 
line: 



Co: 






(12) 



Inspection of Eq. (12) reveals that the numerator repre- 
sents the total ^^^Rn activity passing through the sam- 
pling line from time t^ to t^, while the denominator is the 
corresponding total volume. 

The numerical differences in the mean Co calculated 
by the three calculational approaches were in all cases 
negligible compared to the estimated relative standard 
deviations of the mean (s^) for Co calculated by any one 
of the approaches. These differences, depending on the 
specific standardized addition in question, ranged from 
less than 0.01 % to about 2.0 % to 4.0 % in worst cases. 
The specific case-to-case variation in the precision in 
mean Co may be seen in Tables 2, 3, and 4, which 



Table 1. MIST standardized sample additions 



No. 



Approximate time 
Date— 1991 Hours (GMT) 



Nominal range for the Rn activity concentration 
(Bq ■ m-') 



1 


9 


October 


1900 - 


2300 


2 


10 


October 


0900 - 


1300 


3 


10 


October 


1400 - 


1800 


4 


11 


October 


0900 - 


1300 


5 


11 


October 


1400 - 


1800 


6 


11 


October 


1800 - 


2200 


7 


11-12 


October 


2200 - 


0200 


8 


12 


October 


0200 - 


0600 


9 


12 


October 


1200 - 


1600 





12 


October 


1600 - 


2000 


.1 


12 


October 


2000 - 


0000 


2 


13 


October 


0000 - 


0400 


.3 


13 


October 


1200 - 


1500 


A 


13 


October 


1500 - 


1900 


.5 


13 


October 


1900 - 


2300 



33 to 37 

28 to 30 

29 to 31 
discarded 

5.3 to 6.2 
20 to 23 
3.9 to 5.0 
3.8 to 4.1 
7.2 to 7.9 

12 to 14 
32 to 36 

30 to 33 
15 to 16 

5.4 to 5.8 
2.4 to 2.6 
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Table 2. NIST standard additions for Lab E 





Flow rate /, 




^^^Rn Cone. 
Cs 


Dilution faetor 
D 




"""Rn Cone. 

Co" 


No. 


Mean 


a 


Mean 


a 






(L ■ min"') 


(%) 


(Bq ■ L-') 


(X 10-') 


(%) 


(Bq ■ m-') 


1 


0.335 


0.02 


41.92 


0.854 


0.19 


35.81 




0.333 


0.02 


42.15 


0.796 


0.03 


33.55 




0.332 


0.01 


42.25 


0.778 


0.02 


32.87 




0.328 


0.01 


42.83 


0.865 


0.04 


37.07 


2 


1.200 


0.01 


9.968 


2.919 


0.18 


29.09 




1.201 


0.03 


9.960 


2.804 


0.04 


27.92 




1.194 


0.01 


10.03 


2.874 


0.14 


27.82 




1.186 


0.01 


10.11 


3.014 


0.03 


30.48 


3 


1.173 


0.02 


10.25 


2.909 


0.22 


29.81 




1.175 


0.02 


10.24 


2.792 


0.04 


28.59 




1.166 


0.02 


10.33 


2.946 


0.25 


30.44 




1.175 


0.11 


10.23 


3.052 


0.31 


31.23 


5 


0.985 


0.04 


12.67 


0.422 


0.04 


5.349 




0.981 


0.01 


12.72 


0.430 


0.05 


5.467 




0.979 


0.01 


12.75 


0.454 


0.05 


5.787 




0.990 


0.03 


12.58 


0.460 


0.08 


5.782 


6 


0.768 


0.06 


16.91 


1.298 


0.41 


21.94 




0.770 


0.02 


16.87 


1.211 


0.04 


20.43 




0.773 


0.02 


16.79 


1.294 


0.27 


21.74 




0.772 


0.03 


16.82 


1.300 


0.04 


21.87 


7 


0.645 


0.41 


20.60 


0.243 


3.0 


5.013 




0.664 


0.02 


19.96 


0.207 


1.1 


4.141 




0.668 


0.01 


19.82 


0.199 


0.06 


3.946 




0.670 


0.01 


19.73 


0.210 


0.10 


4.143 


8 


0.676 


0.01 


19.54 


0.208 


0.04 


4.064 




0.669 


0.02 


19.78 


0.201 


0.14 


3.978 




0.671 


0.01 


19.71 


0.195 


0.04 


3.833 




0.671 


0.01 


19.70 


0.198 


0.15 


3.898 


9 


1.051 


0.12 


11.72 


0.627 


0.63 


7.350 




1.042 


0.02 


11.84 


0.607 


0.03 


7.181 




1.019 


0.05 


12.16 


0.619 


0.19 


7.524 




1.026 


0.02 


12.07 


0.654 


0.07 


7.885 


10 


0.958 


0.03 


13.08 


1.025 


0.15 


13.41 




0.944 


0.01 


13.32 


0.930 


0.03 


12.39 




0.944 


0.01 


13.31 


0.952 


0.05 


12.67 




0.953 


0.02 


13.17 


1.041 


0.11 


13.70 


11 


0.710 


0.02 


18.49 


1.961 


0.02 


36.25 




0.710 


0.01 


18.50 


1.863 


0.26 


34.46 




0.713 


0.02 


18.40 


1.764 


0.06 


32.47 




0.717 


0.01 


18.29 


1.831 


0.25 


33.49 


12 


0.905 


0.01 


14.00 


2.391 


0.02 


33.48 




0.909 


0.01 


13.92 


2.278 


0.12 


31.71 




0.916 


0.02 


13.80 


2.051 


0.03 


28.30 




0.907 


0.01 


13.95 


2.124 


0.03 


29.64 


13 


0.916 


0.10 


13.80 


1.072 


0.16 


14.79 




0.920 


0.02 


13.73 


1.142 


0.14 


15.68 




0.922 


0.01 


13.69 


1.144 


0.02 


15.66 


14 


1.073 


0.15 


11.42 


0.505 


1.33 


5.770 




1.078 


0.01 


11.37 


0.477 


0.03 


5.422 




1.079 


0.01 


11.36 


0.499 


0.21 


5.671 




1.081 


0.01 


11.33 


0.513 


0.03 


5.811 


15 


1.146 


0.03 


10.56 


0.242 


0.58 


2.556 




1.145 


0.01 


10.56 


0.225 


0.06 


2.371 




1.143 


0.01 


10.59 


0.225 


0.85 


2.383 




1.138 


0.01 


10.65 


0.232 


0.10 


2.465 



' Relative standard deviation of the mean. 

' The relative "overall uncertainties" ranged from 6.0 % to 13.0 % (see Table 7). 
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Table 3. NIST standard additions for Lab A 





Flow rate/4 




^^^Rn Cone. 
Cs 


Dilution factor 
D 




"""Rn Cone. 

Co" 


No. 


Mean 


a 


Mean 


a 






(L • min"') 


(%) 


(Bq • L-') 


(X 10-') 


(%) 


(Bq • m-') 


1 


0.332 


0.02 


42.30 


0.825 


0.21 


34.91 


2 


1.195 


0.02 


10.02 


2.915 


0.20 


29.20 


3 


1.172 


0.03 


10.26 


2.935 


0.21 


30.11 


5 


0.984 


0.02 


12.68 


0.480 


1.9 


6.09 


6 


0.771 


0.02 


16.84 


1.338 


2.0 


22.53 


7 


0.662 


0.11 


20.02 


0.217 


1.0 


4.35 


8 


0.672 


0.01 


19.68 


0.206 


1.2 


4.06 


9 


1.035 


0.05 


11.94 


0.643 


0.89 


7.68 


10 


0.950 


0.02 


13.22 


1.003 


0.57 


13.27 


11 


0.713 


0.01 


18.42 


1.900 


0.88 


35.01 


12 


0.909 


0.01 


13.92 


2.232 


0.33 


31.06 


13 


0.919 


0.03 


13.74 


1.149 


1.1 


15.78 


14 


1.078 


0.04 


11.37 


0.510 


0.86 


5.80 


15 


1.143 


0.01 


10.59 


0.238 


1.9 


2.52 



Relative standard deviation of the mean. 
' The relative "overall uncertainties" ranged from 6.0 % to 13.0 % (see Table 7). 



Table 4. NIST standard additions for Lab D 





Flow rate/4 




^''Rn Cone. 


Dilution factor 
D 




'"'Rn Cone. 
r ^ 


No. 


Mean 


a 


Mean 


a 






(L • min"') 


(%) 


(Bq • L-') 


(X 10-') 


(%) 


(Bq • m-') 


1 


0.332 


0.02 


42.21 


0.792 


0.27 


33.43 


2 


1.198 


0.02 


9.99 


2.844 


0.26 


28.41 


3 


1.174 


0.02 


10.25 


2.823 


0.26 


28.93 


5 


0.983 


0.02 


12.70 


0.490 


3.2 


6.22 


6 


0.771 


0.02 


16.85 


1.291 


3.0 


21.76 


7 


0.666 


0.01 


19.88 


0.205 


0.81 


4.08 


8 


0.671 


0.01 


19.72 


0.202 


2.3 


3.98 


9 


1.040 


0.04 


11.88 


0.619 


0.94 


7.35 


10 


0.944 


0.01 


13.32 


0.964 


0.94 


12.84 


11 


0.714 


0.02 


18.39 


1.814 


1.09 


33.36 


12 


0.912 


0.02 


13.88 


2.104 


0.31 


29.20 


13 


0.917 


0.09 


13.78 


(0.612)" 


1.2 


(8.43)" 


14 


1.079 


0.01 


11.36 


0.488 


0.91 


5.54 


15 


1.143 


0.01 


10.58 


0.234 


3.6 


2.48 



' Relative standard deviation of the mean. 

' The relative "overall uncertainties" ranged from 6.0 % to 13.0 % (see Table 7). 

' This standard addition (#13) for Lab D was normalized by the ratio 0.5583 to account for a partial standard addition during the measurement 
interval. The standard addition was made for only the last 67 min of the 120 min measurement interval. 
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provide the calculated results of the standardized addi- 
tions for the three participating laboratories. Estimates 
of Sm for Co may be obtained directly from Sm for the 
mean dilution factor D since the latter is the singularly 
dominant source of variability in Co- These estimates of 
Sn range from about 0.03 % to 3.0 % for Lab E, and 
from about 0.2 % to 3.6 % for Lab A and Lab D. The 
wide variations merely reflect the actual wide variations 
in the flow rates (mainly /i) from case to case as will be 
discussed shortly. Nevertheless, any differences due to 
the chosen calculational approach for a given case were 
clearly reflected and embodied within the estimated Sm 
for that case. 

The final calculational approach chosen for the inter- 
comparison of standardized additions consisted of a 
combination of the first two approaches outlined above. 
It must be emphasized, however, that this choice was 
somewhat arbitrary since all of the approaches gave 
essentially identical results. None of the results and con- 
clusions of the intercomparison would change as a result 
of a different calculational choice. For the adopted 
method chosen, the mean Cs was obtained from a mean 
/4 [Eq. (10)] as in the first approach; the mean dilution 
factor D [Eq. (7)] was derived as in the second 
approach; and the mean Co was just the product of these 
values of Cs and D. These are the values tabulated in 
Tables 2, 3, and 4. 

As indicated in the tables, the ^^^Rn activity concen- 
trations Co in the standardized sample additions ranged 
from approximately 2.5 Bq • m"^ to 35 Bq • m"^. 

Perusal of the actual flow-rate measurement data for 
a typical standardized sample addition may be helpful in 
understanding the experimental aspects of the intercom- 
parison and the quality of the results. Figures 4 through 
7 give the flow -rate measurement data for/4,/3,/5, and 
/i, respectively, obtained during standardized addition 
#11. The illustrated data consists of nearly 1 440 simul- 
taneous measurements of the four flow rates, taken 
approximately every 10 s over an interval of 4 h. 

The means for the relatively stable flow rates /4,/3, 
and /s were somewhat over determined and are very 
precise. Tables 2, 3, and 4 give the estimated relative 
standard deviations of the mean {Sm) for/4, and these 
values of Sm are typical of those for/3 and/5 as well. In 
fact, even these values are somewhat misleading for the 
inherent flow-rate variability in the middle of a standard 
addition. They are sometimes strongly influenced by 
abrupt changes in the flow rate that were made over the 
first 0.5 min to 2 min in adjusting the dilution factor in 
going from one standard addition to the next. Obviously, 
having a time delay between sequential additions would 
have obviated this shortcoming, but the rather costly 
nature of this intercomparison imposed time constraints. 
This kind of adjustment between standard additions is 



clearly shown in Fig. 5 where flow rate/3 was adjusted 
to vary the dilution factor from standard addition # 10 to 
#11. Nevertheless, the data of Figs. 4 through 7 nicely 
illustrate the typical variations in flow over the course of 
a standard addition. The flow-rate data also often exhib- 
ited some slight, but gradual systematic drifts in flow 
rate over time. This effect is quite apparent in the data of 
/3 and/5 shown in Figs. 5 and 6. This is not believed to 
be a result of any electrical or signal "drift" in the 
flow meters themselves at constant flow, but rather is 
believed to have arisen from gradual changes in the 
settings of the flow control valves in the manifold due to 
mechanical vibrations at the site. 

In contrast to the flow rates /4,/3, and/5, the data for 
flow rate/i, shown in Fig. 7, exhibit pronounced shifts. 
These marked changes arise as a result of the on and off 
cycling of the sampling pumps that were part of the 
measurement instruments for the three participating lab- 
oratories. The flow rate/i at any time was the sum of the 
flow rates from the three laboratories' pumps. Mean 
flow rates/] calculated for each individual standardized 
addition over the course of all the additions ranged from 
approximately 365 L • min"' to 450 L • min"', with 
estimated relative standard deviations of the mean {Sm) 
of essentially the same magnitude as that for the dilution 
factors D (see Tables 2, 3, and 4). The variability in/i 
was the predominant contribution to the variability in 
the dilution factor D . The Lab E pump, which rapidly 
cycled on and off (over approximately 40 s) every hour 
was dominant, contributing about 365 L • min"' to 
395 L • min"' of the total flow rate. The Lab D pump 
with a flow rate of about 40 L • min"' to 45 L • min"' 
sampled for 2 h, but at irregular intervals of about every 
4 h to 5 h. The Lab A pump ran continuously with a 
flow rate close to 40 L • min"' that increased to about 
44 L • min"' for 20 s every 10 min. The 4 intervals 
labelled Ella through El Id in Fig. 7 mark the pump 
sampling intervals for Lab E. The rapid drops in the 
flow rate as the Lab E pump shut down and then 
returned is very pronounced in the data, but is barely 
perceptible in Fig. 7 since each cycle only affects 3 or 
4 out of the 1440 data values shown. The interval 
marked Dll in Fig. 7 represents the sampling interval 
for Lab D. 

The dilution factor D at each 10 s interval for the 
standard addition # 1 1 flow data (Figs. 4 through 7) is 
given in Fig. 8. Of necessity, it exhibits the same marked 
discontinuities as the/i data. It is apparent, for example, 
that the dilution factor D (and hence Co) is quite stable 
for the Lab Ella and El Ic intervals; and is somewhat 
less so for the Lab D 1 1 interval because of the two 
disruptions by the Lab E pump. The intervals El lb and 
El Id exhibit pronounced step functions in the dilution 
factor as the Lab D pump went on and off. These evident 
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Fig. 4. Typical flow rate/4 data as measured by meter F4 approximately every 10 s over an interval of 4 h for 
standard addition #11. 
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Fig. 5. Typical flow rate/3 data as measured by meter F^ approximately every 10 s over an interval of 4 h for 
standard addition #11. The discontinuity in flow at the outset (from about 0.73 L • min"' to 0.67 L • min"') was 
an adjustment in the flow rate to change the dilution factor from standard addition #10 to #11 (see text). 
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Fig. 6. Typical flow rate/5 data as measured by meter F5 approximately every 10 s over an interval of 4 h for 
standard addition #11. 
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Fig. 7. Typical flow rate/i data as measured by meter Fi approximately every 10 s over an interval of 4 h for 
standard addition #11. The discontinuities are discussed in the text. 



Q 

O 

u 

a 



e 
o 



3 

a 



0.0025 
0.0023 - 
0.0021 - 
0.0019 - 
0.001 7 



D11 



Ella 



E11b 



Ellc 



Elld 



1 — I — r 



jSSf!i»i>*gfPOI 



0.0015 
285.82 



.J ! L. 






-I L 



' j , a fc' » '''*M ' '» ' ** ef 



J I L 



285.86 285.90 285.94 

Time GMT (1991 Julian day) 



285.98 



286.02 



Fig. 8. Dilution factor D as a function of time over the interval of 4 h for standard addition #11. The data values 
were calculated every 10 s using the flow rate data of Figs. 4 through 7. 
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and very different kinds of variations in the dilution 
factor are also manifested by the magnitude of their 
respective estimated relative standard deviations of the 
mean (s^)- 

Consider the respective cases for the three labora- 
tories. 

Table 2 indicates that for the stable Ella and Ellc 
intervals, Sm for dilution factor D is 0.02 % and 0.06 %; 
while for intervals El lb and El Id, with the step 
function in D, Sm is 0.26 % and 0.25 %, respectively. 
Intuitively, one might expect s^ for mean D for these 
cases with a drastic step function in D to be even larger 
than that given. This estimated Sm, however, is easily 
derived and verified from its component parts. For a 
step function with mean OTi and Smi at the low plateau 
having a time duration At,, and with mean m2 and s^i at 
the high plateau for time duration At2, the mean m over 
the entire (Ati + At2) duration is approximately 



Ati(mi) + At2(m2) 
(Ati + Ah) 



(13) 



and the relative standard deviation of the mean for m , 
obtained by appropriately adding the variances of rrii 
and m2, is 



The mean dilution factor D for interval Dl 1 has a s^ 
for £) of 1 . 1 % (Table 4). Although it is not as apparent 
in the illustrated data of Fig. 8, its magnitude is strongly 
influenced by the two abrupt flow discontinuities when 
the Lab E pump went off and on. 

Over the entire 4 h interval of Fig. 8, as applicable for 
calculation of the Lab A dilution factor mean D , Sm for 
D was 0.9 % (Table 3). 



5. Confirmatory Measurements 

The dimensionless dilution factors D [Eq. (7)], used 
to obtain the values of Co for the intercomparison, 
consist of flow rate ratios, and therefore are more depen- 
dent on the relative flow meter responses rather than on 
individual flow meter calibrations per se . Although the 
matched flow meters F4, F3, and F5 in the NIST manifold 
were previously intercompared and also presumably 
based on similar and relatable calibrations by the 
manufacturer (Matheson), the design of the manifold 
(Fig. 1) and the simultaneous flow-rate measurements 
allowed a direct, in situ intercomparison of the flow 
meters at the test site. 

Inspection of Figs. 1 and 2 indicates that when/c = 0, 
then, in principle, (f^ +f^) =fs. This provided the oppor- 
tunity to intercompare the flow-meter responses for the 
following conditions 



A?, 



(At,+At2) 




A?, 



(Ati+At2) 




(14) 






h 



/s 



with 



with 



/3 = 



/4 = 



(15) 



Clearly, the mean m is proportionately-distant between 
the low and high plateaus of the step function, and the 
negative residuals at the low plateau are counterbalanced 
by all of the positive residuals at the high plateau. 
Furthermore, the step function in dilution factor D is 
inherently bimodally distributed (i.e., a bivariate distri- 
bution) in these cases, and s^ for the overall mean m is 
an indescriptive and non-robust statistical estimator of 
the variability in D . The estimated Sm, nevertheless, if 
anything, overestimates the true dispersion in the two 
bimodal means. In view of the invariance (within 
statistical variations) of the three different calculational 
approaches to obtain a mean Co (Sec. 4), it is clear that 
even with the step function in D (and thereby in Co), the 
adopted treatment is adequate. 



/4+/3 



All three conditions were tested in Bermuda during the 
course of the intercomparison. The results are given in 
Table 5. Mean values for the flow rates were obtained 
from the 10 s measurement data averaged over periods 
of from 20 min to nearly 60 min. All of the estimated 
relative standard deviations of the mean for the flow 
rates were in the range 0.01 % to 0.04 %. As indicated, 
all of the flow-ratio intercomparisons are within ± 1 %, 
which is well within the range of the expected statistical 
variations for meters with an assumed relative standard 
uncertainty of 0.5 % (see Sec. 3 discussion). 
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Regrettably, it was not possible to independently 
verify by intercomparison the relative response of flow 
meter Fi to the others. The repercussions of this short- 
coming as it affects the results and conclusions of 
the intercomparison of standardized additions will be 
addressed in Part II of this series [1]. 

The conformity of the flow rates at sequential 
locations in the streamlines of the manifold also demon- 
strated the absence of any serious leaks in the manifold's 
construction. 

Two types of in situ confirmatory tests of the opera- 
tion of the NIST manifold to deliver known ^^^Rn activ- 
ity concentrations at the site in Bermuda were per- 
formed. Both involved collecting "grab samples" from 
the manifold using glass sampling bulbs identical to 
those used to perform the ^^''Ra source calibration. The 
samples were returned to NIST for assay, again using a 



PIC measurement procedure identical to that used in 
calibrating the radium source (see Sec. 3). The assayed 
^^^Rn activity concentrations in the samples were decay 
corrected to the time of sample collection, and this 
corrected concentration Ct was compared to that 
concentration Cp predicted from the operation of the 
manifold. 

The first kind of test involved samples collected at the 
sampling port directly downstream from the radium 
source (see Fig. 1). The activity concentrations in these 
samples are expected to be that given by the source 
calibration [Cs of Eq. (10)] as a function of the flow rate 
/4. Three such samples were collected during standard 
additions #7, #10, and #14, with the results summa- 
rized in Table 6. An additional two samples were 
collected, but were unfortunately destroyed in transit to 
NIST. 



Table 5. In situ intercomparison of the flow meters F4, F3, and F5 



Test 
conditions 



/3 

Mean 



Flow rates (L ■ min ') 



(%) 



/4 

Mean 



(%) 



(Mean) 



(%) 



Flow 
ratio 



/4 = 
/4=/5 



/3+/4=/5 



0.770 



0.681 



0.02 



0.03 







0.766 


0.03 


- = 1.005 
h 


0.793 


0.01 


0.797 


0.03 


/4 

- = 0.995 

h 


0.919 


0.01 


0.927 


0.01 


= 0.991 


1.078 


0.04 


1.076 


0.02 


= 1.002 


0.713 


0.01 


1.400 


0.02 


/3./4__^ 



0.996 



' Estimated relative standard deviation of the mean. 



Table 6. Results of the confirmatory measurements for Rn activity concentration collected in situ at Bermuda and assayed at NIST 











i- -. 












/4 


Cs^ 


/3 


/4 


Cp" 




Ct^ 




Sample 


(L ■ min"') 


(Bq ■ L-') 


(L ■ min"') 


[/4+/3J 


(Bq ■ L-') 




(Bq ■ L-') 


Cp/Ct 


NIO 


0.9572 


13.10 


(0) 


(1) 


13.10 




13.27 


0.9873 


N14 


1.0840 


11.29 


(0) 


(1) 


11.29 




11.43 


0.9878 


N7(b) 


0.6924 


19.03 


(0) 


(1) 


19.03 




19.43 


0.9792 


ND 


0.4715 


29.06 


0.6948 


0.4043 


11.75 




11.69 


1.0049 


NF 


0.3932 


35.32 


0.6983 


0.3602 


12.72 




13.01 


0.9779 












Number 




5 












Cp/Ct 


mean 


0.9874 












Sm 


(%) 


d 


0.49 



" See Eq. (10). 

Cp is the predicted ^^^Rn activity concentration given by Cp 

Ct is the assayed, decay-corrected ' 
refer to text for details. 
Relative standard deviation of mean in percent. 



l/4/(f4+/3)]Cs. 

Rn activity concentration at the "grab sample" time of collection, corresponding to that expected for Cp; 
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The second kind of test was performed at the con- 
clusion of the standardized additions when the NIST 
manifold was disconnected from the main sampling 
line. Two samples were collected directly at the outlet of 
the manifold (see Fig. 1) after meter F5 and just as it 
would have entered the main line before meter Fi. This 
confirmatory test was made to demonstrate the efficacy 
of mixing the/4 and/3 flows, and the (/4 +/3) dilution 
calculation. In this case, the predicted concentration Cp 
is given by Eq. (3), where Cs is diluted by the factor 
f^/if^ +f?)- The results for theses two samples, labelled 
ND and NF, are also presented in Table 6. 

The mean flow rates for fa, and f^ used to obtain 
the predicted concentration Cs and diluting factor 
fJifA+f^) were obtained from the 10 s flow -rate data 
averaged over intervals of approximately 10 min that 
immediately preceded the closing of the sample bulbs. 
The estimated relative standard deviation of the mean 
for/4 and/3 was, in all cases, less than 0.01 %. The mean 
values averaged over shorter intervals (down to 3 min) 
and longer intervals (up to 20 min) were invariant. 

Comparison of the assayed concentrations Ct and 
predicted values Cp in Table 6 indicates an overall agree- 
ment to slightly over 1 %, with no significant differ- 
ences between the two kinds of samples. The mean 
value of Cp/Cx = 0.987 has an estimated relative 
standard deviation of the mean (^m) of 0.49 %. This 
uncertainty is roughly the same magnitude as that 
expected for the derived uncertainty in Cs as given ear- 
lier (Sec. 3). As indicated, a very slight negative bias in 
Cp/Cx may be suggested, particularly for the three 
direct Cs samples, "iet, even if this bias exists, the statis- 
tical variations and small number of samples does not 
allow its confirmation. It should be noted that the Cp/Cx 
comparison given in Table 6 neglected any contribution 
of radon activity from ambient air. In actuality, a com- 
parison of (Cp + Ca)/Cx (where Ca is the ambient ^^^Rn 
activity concentration at the time of sample collection) 
would be more valid. However, Ca is believed to be 
^ 0.003 Bq ■ L"' for any of the samples, and should 
therefore have a negligible effect (^ 0.03 % of Cx) on 
the results. This limit is based on the maximum ambient 
concentration that was measured by any of the partici- 
pating laboratories at any time during the intercompari- 
son (see Ref. [1]). 



6. Uncertainty Analysis 

A complete analysis of the measurement uncertainties 
for the ^^^Ra source calibration and ^^^Rn activity 
concentrations in the standardized sample additions is 
outlined in extenso in Table 7. 



The uncertainty analysis procedure follows the 
normal conventions of the NIST Radioactivity Group 
which for the most part are compatible with those 
adapted by the principal international metrology 
standardization bodies [4,5]. All individual uncertainty 
components are expressed in terms of estimated (exper- 
imental) standard deviations (or standard deviations of 
the mean where appropriate) or quantities assumed to 
correspond to standard deviations, irrespective of the 
method used to evaluate their magnitude. A propagated 
or ' 'combined standard uncertainty' ' is expressed as an 
estimated standard deviation which is equal to the 
positive square root of the total variance obtained by 
summing all variance and covariance components, 
however evaluated, using the law of propagation of un- 
certainty for the specific mathematical function given 
by the model of the measurement procedure. By 
convention in this laboratory (at the time of this inter- 
comparison),^ the combined standard uncertainty is 
expanded by a "coverage factor" of 3 to obtain an 
"expanded uncertainty" (or "overall uncertainty" 
[sic]) which is assumed to provide an uncertainty inter- 
val having a high level of confidence of roughly 95 % to 
99%. 

The analysis of the propagated uncertainty in the con- 
centration Cs exiting the radium source as a function of 
the flow rate/4 was discussed previously in Sec. 3. 

To derive the uncertainty in the concentration Co for 
the standardized sample additions from the uncertainty 
in Cs, one must consider the additional uncertainty 
components due to: (1) the mean for flow rate/4 for the 
given addition; (2) the mean dilution factor D ; (3) the 
flow-meter calibrations; and (4) the timing of the addi- 
tion which determines the start and stop points of the 
numerical calculation of the fa, and D means. 

The statistical estimators (s^) for f^ and D were 
treated and discussed at considerable length in the 
example for addition #1 1 in Sec. 4. The values listed in 
Table 7 have wide ranges exceeding an order of magni- 
tude. The reasons for the wide variations were also 
addressed in Sec. 4. 

The uncertainties in the flow-meter calibrations, 
again expressed as relative standard deviations, were 
obtained by taking one half of the manufacturers' stated 
calibration accuracies since the latter are presumed to 
correspond to a high level of confidence. The uncer- 
tainty in the timing of the standardized addition is also 
highly variable depending on the specific case. The 
uncertainties {A) in the actual start {t^ and stop {t^ 
times per se are not critical, but rather it is the uncer- 
tainty in calculating the means /4 and D over the interval 



NIST has since uniformly adopted a policy of using a coverage factor 
of yt = 2. 
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Table 7. Analysis of uncertainties for the ^^^Ra source calibration and ^^^Rn concentrations in the standardized additions 



Source of uncertainty 

Precision of radon measurement in 
sample bulbs by PIC (5,) 

PIC calibration (5,) 

Volume sample bulbs (5,) 

Radium and radon decay corrections (5,) 



Radon concentration in sample \ 5s 

Precision flow rate during sampling ( 5,- ) 
Flow meter calibration (5,) 



2,- 5? 



Flow rate 



S, sf 



Fit (regression) of radon concentration as a function 
of flow rate (5r) 

Radon concentration (Cs) from source at given flow 



rate /i 



Ss + Sf + Sr 



Timing errors for standardized addition (5*) 

Overall precision in diulution factor from flow rate 
measurement (St) 

Calibration flow meters Fj, F4, and F5 (5*) 

Calibration flow meter Fi (5*) 



Dilution factor \ 5d = A / St Sl 
Precision flow rate/4 for std. addition (5/^) 
Radon concentration (Co) in std. addition 



(excluding ambient cone.) dc„ = A / 5^ + Sj + 5}^ 

Radon concentration (Ca) for ambient air (Scf) 

Radon concentration (Ci) for std. addition 
(including ambient cone.) 




g) «^..[g) ^h 



For Co / Ca = 100 
For CqI Ci^= 25 
ForCo/CA= 10 



Component 


Propagated 


Typical "overall 


uncertainty " 


uncertainty " 


uncertainty' ' 


0.05 to 0.2 






Q.r 


— 


— 


0.2 to 0.5 


— 


— 


<0.01 


— 


— 




0.73 to 0.88 


2.2 to 2.6 



0.02 
0.5'' 



0.77 to 2.4 



0.1 to 0.3 

0.2 to 3.0 
0.5'' 
1.4" 



0.1 to 0.4 



Assumed 50 



0.5 



1.2 to 2.6 



1.6 to 3.5 



2.0 to 4.4 



2.0 to 4.4 
2.7 to 4.6 
4.9 to 6.1 



3.5 to 7.i 



6.0 to 13 



6.1 to 13 

8.2 to 14 
15 to 18 



The estimated component and propagated uncertainties are assumed to correspond to an approximate relative standard deviation (or standard 
deviation of the mean) expressed in percent. 

The "overall uncertainty" (or expanded uncertainty with coverage factor i = 3) is taken to be three times the propagated uncertainty in percent. 
" See Refs. [2] and [3]. 



Corresponds to one half the vendor's (Matheson) stated calibration accuracy. 
' Corresponds to one half the vendor's (TSI, Inc.) stated calibration "tolerance.' 
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from (fa ± A, J to (f^ ± A, J that is of importance. The 
additions for Lab D were largely free of this uncertainty. 
Their measurement intervals were (with one exception) 
completely overlapped by the 4 h or 3 h standard addi- 
tion intervals. Shifting fa and f^ by even several minutes 
in either direction has almost negligible effects 
(< 0.1 %) on the results of the numerical integrations for 
the Lab D additions. The Lab A additions utilized the 
calculational results for the entire 4 h or 3 h sample 
intervals. It is conceivable that the clock times for their 
sampling intervals of 1 h differed from the assumed 
fa and fz times by 1 min to 2 min. This presented the real 
possibility of having overflows of radon from one 
presumed standardized addition interval to another in 
the case of adjacent addition intervals. Perhaps, this kind 
of source of inaccuracy should be considered to fall 
within the category of a ' 'blunder' ' rather than an eval- 
uatable uncertainty component, "iet, considering the 
quality and nature of their "smoothed" measurement 
results and the somewhat equivocable, attendant analysis 
of their data that was thereby required [1], this uncer- 
tainty component was truly negligible in terms of any 
interpretations of the measurement results and conclu- 
sions of the intercomparison. The abrupt flow interrup- 
tions by the pump for Lab E served as an internal timer 
to determine the start and stop times of the additions 
(see Sees. 2 and 4) so, in principle, their additions were 
nearly perfectly synchronized in terms of real clock 
times. However, even the small 15 s to 20 s A,^ and A,^ 
intervals, which had accompanying large fluctuations in 
the flow rates just at the fa and f^ times, resulted in the 
largest uncertainties in f^ and D. These approached 
0.3 % of /t in worst cases. 

The above component uncertainties may be propa- 
gated, as done in Table 7, to form the uncertainty in Co. 
It ranges from about 2.0 % to 4.4 % of Co, or about 6 % 
to 13 % at the "overall uncertainty" level. This is the 
uncertainty associated with the ^^^Rn activity concen- 
tration in the standardized sample additions, exclusive of 
any contributions from ^^^Rn in the ambient air. Consid- 
ering the uncertainty associated with the uncertainty 
Scq itself (as well as the largely unknown degrees of 
freedom and the unknown probability distributions for 
the component quantities), individual values of Sc„ were 
not propagated for each of the 83 standardized sample 
additions listed in Tables 2, 3, and 4. That kind of arith- 
metic exercise was deemed to have little significant 
value. The detailed uncertainty analysis presented here 
resulted in a range of 6 % to 13 % for the "overall 
uncertainty" in Co, which justifiably may be considered 
to be in the general range of approximately 10 %. 

Now, to treat the uncertainty in the ^^^Rn activity 
concentration C\ in the main sampling line as sampled 
by the participating laboratories, one must consider the 



contribution due to the ambient concentration, Ca- 
Recalling that Ci = Co + Ca [Eq. (5)], the uncertainty in 
C\ is given by the sum of the variances Uc^ = Uc„ + Uc^, 
which yields 




3c„ 



(16) 



when the uncertainties (S) are expressed as relative 
standard deviations. Obviously, the uncertainty in 5c, 
can vary extensively depending on the magnitude of the 
concentration ratio CJCa and the uncertainty in the 
ambient concentration (ScjJ- Table 7 illustrates this for 
the hypothetical case of an assumed ambient concentra- 
tion uncertainty (Sc,) of 50 % at three CJCa ratios. For 
Co/Ca= 10, the uncertainty 5c, increases by about a 
factor of 1.5 to 2.5 times 8co, while for CJCa = 100, 
there is negligible difference between 5c, and 5c „• 
Throughout the course of the intercomparison for all 
standardized sample additions, the ratio CJCa was 
estimated to lie in the range 10 < CJCa < 1000 [1]. 



7. Conclusions 

This work provided a standardized reference basis for 
the first international intercomparison of ^^^Rn detectors 
used in marine-atmospheric studies. The work went 
beyond serving the needs of this particular intercom- 
parison. More importantly, it also demonstrated the 
broader utility of the calibration protocol and the 
methodology for the standardized sample additions that 
were developed for it. 

The intercomparison not only provided a basis for 
comparing the measurement results and performance of 
different instruments of various participating laborato- 
ries, but more so, provided a common reference to these 
laboratories and provided the possibility for an in situ 
intercalibration that could be related to U.S. national 
^^''Ra standards. Most environmental measurement 
intercomparisons of field instruments in actual use 
merely rely on evaluating the relative performance of 
the participants, or some comparison to the pooled 
results. This exercise demonstrated, for the very first 
time, the capability of providing a standardized refer- 
ence basis even for such low-level, field-measurement 
intercomparisons. 

The standardized sample additions of known, 
but undisclosed ("blind") ^^^Rn activity concentration 
used for this intercomparison ranged from approxi- 
mately 2.5 Bq • m"^ to 35 Bq • m"^, and had "overall 
uncertainties" (that can be related to national standards) 
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in the general range of approximately 10 % at an 
assumed three standard deviation uncertainty interval. 
As will be presented and discussed in the second paper 
of this series [1], the ^^^Rn activity concentrations in 
Bermudian ambient air over the course of the intercom- 
parison are understood to range from < 0.01 Bq • m"^ to 
roughly 2 Bq • m"^. Thus, there was nearly a complete 
and compatible overlap between the concentrations in 
the standardized sample additions and in ambient air. 

The developed methodologies presented here could, 
of course, be adopted with slight modifications to cover 
other ^^^Rn concentration ranges and other applications, 
and could be employed in many other types of ^^^Rn 
environmental, field-measurement intercomparisons. 
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